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ABSTRACT 

The  theory  of  shockwave  focusing  (developed  by 

Silbiger^)  is  employed  in  an  attempt  to  predict  the  extent 

of  focusing  of  explosive  pulses  at  caustics  as  observed 

2 

experimentally  by  Barash. 

Assuming  a  viscous  medium  end  a  physically  reasonable 
model  for  the  explosive  pulse,  an  order  of  magnitude  com¬ 
parison  with  experiment  is  obtained.  The  omission  of 
viscosity  from  the  mathematical  model  leads  to  unaccept¬ 
ably  large  amplification  at  the  caustic. 
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i.  introduction 


Because  the  speed  of  sound  varies  from  point-to-point  in  most  bodies  of 
water,  shockwaves  propagating  under  water  do  not  travel  in  straight  lines 
but  rather  in  curved  paths.  Under  certain  conditions  and  in  certain  regions 
the  paths  converge,  or  focus,  and  sound  energy  which  had  previously  been  dis¬ 
tributed  over  a  large  volume  of  water  is  concentrated  into  a  small  volume. 
Regions  in  which  this  effect  occurs  are  called  convergence  zones.  Conver¬ 
gence  zones  consist  of  a  point,  line  or  surface  -  called  the  caustic  -  where 
the  focusing  is  maximum,  surrounded  by  a  narrow  region  in  which  the  focusing 
effect  diminishes  with  distance  from  the  caustic. 

For  submerged  sources  of  explosive  sound,  the  pressure  pulses  appearing 
at  the  caustics  display  marked  amplification  of  the  peak  value  and  distor¬ 
tion  of  the  waveform. 

A  theory  to  describe  the  propagation  ana  focusing  of  transient  pulses  in 
convergence  zones  has  been  proposed  by  Silbiger.1  The  theory  is  based  on 
the  fundamental  assumptions  of  the  validity  of  linear  acoustics  and  the.  ab¬ 
sence  of  viscosity  in  the  propagating  medium.  Modifications  to  incorporate 
viscous  effects  into  the  theory  can  be  made  and  are  discussed  in  a  later  part 
of  this  report.  Non-linear  effects  can  also  be  included  in  the  theory;  how¬ 
ever,  they  are  not  considered  in  this  report. 

By  making  a  series  of  mathematical  approximations,  the  theory  can  be  shown 
to  reduce  to  the  results  of  geometric  acoustics,  in  regions  off  caustics  and 
not  in  shadow  zones.  Through  a  modification  of  geometric  acoustics,  the 
theory  yields  predictions  of  the  pressure  on  the  caustic  itself  and  in  the 
shadow  zone  near  the  caustic.^ 
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The  goal  of  this  report  is  to  report  results  of  a  numerical  comparison 

2 

of  the  Silbiger  theory  with  the  experimental  results  of  Barash. 

The  numerical  comparison  with  experiment  indicates  that  an  order  of  mag¬ 
nitude  prediction  of  the  peak  amplification  of  a  focused  pulse  is  obtained 
assuming  a  viscous  medium  and  zero  rise  time  source  pulse. 

II.  COMPARISON  OF  THEORY  AND  EXPERIMENT 

A.  Experimental  Sound  Speed  Profile  and  Sound  Sources 
2 

Barash  observed  the  focusing  of  underwater  shocxwaves  produced  by  deto¬ 
nating  pentolite  explosives  in  a  quarry  with  a  sound  speed  profile  as  shown 
in  Fig.  1.  In  one  of  his  experiments  the  sound  source  was  placed  at  a 
depth  of  5°  ft  for  which  the  sound  rays  (determined  by  Snell's  law)  are  shown 
in  Fig.  2.  A  caustic  surface  is  evident.  Barash  observed  the  pressure  pulse 
produced  at  the  caustic  by  the  detonation. 

We  shall  be  concerned  with  reproducing  theoretically  the  experimental 
results  for  a  particular  combination  of  experiments]  parameters:  charge 
weight,  0.122  lb;  charge  depth,  50  ft;  and  a  horizontal  ra^ge  (to  the  obser¬ 
vation  point  on  the  caustic)  of  300  ft,  at  which  the  caustic  depth  was 
29.2  ft. 

3.  Mathematical  Model  of  Experiment 
1.  Modeling  the  Source  Pulse 

A  fundamental  assumption  of  the  Silbiger  theory  of  shockwave  focusing 
is  that  non-linear  effects  are  unimportant.  Tnese  can  arise  from  large  am¬ 
plitude  pressure  waves  so  we  chose  to  consider  the  smallest  explosive  weight 
experimentally  employed  ( 0.122  lb)  in  the  expectation  that  linear  acoustics 
theory  will  be  valid. 


-2- 


The  pressure  pulse  produced  in  the  vicinity  of  a  submerged  explosive 

is  usually  said  to  have  a  negligible  rise  time  and  an  exponential  decay,  the 

decay  constant  and  peak  value  being  functions  of  the  charge  weight  and  type, 
3 

and  the  range.  For  pentolite  charges  in  a  non-re fra ctive  medium,  we  have 


p(R,t)  =  2.25x10 


<&■ 
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■(t-R/c) 


>o  (t  <  S)  (1) 

where  w  is  the  charge  weight  in  pounds,  R  is  the  range  in  feet  and  c  is 
the  velocity  of  sound,  t^,  the  decay  constant,  is  given  by 


-  58  w1^  (2L_) 


1/3-0.22 


(micro-secs) 


(2) 


A  "negligible"  rise  time  presumably  means  a  rise  time  much  smaller 
than  experimental  equipment  can  detect.  As  of  yet,  investigators  in  the 
field  of  underwater  detonations  have  not  achieved  a  measurement  of  the  rise 
tine  close  to  the  explosion;  the  generally  quoted  (and  use*)  value  is  zero. 
Clearly  there  is  some  non- zero  yet  very  small  rise  tine;  about  the  only 
quantitative  statement  that  can  oe  made  is  that  the  rise  tine  must  be  much 
shorter  than  an  observable  decay  time. 

In  Siibiger's  notation,  the  pressure  radiated  by  a  transient  source 
into  an  inhomogeneous  medium  obeys  the  following  wave  equation 

y2P(x,y,2,t)  -  ~  p(x,y,z,t)  =  -  ;(x)6(y)6(z)  ?(t)  (3) 

c  (z)  St 
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V 

i, 

C: 


where  close  to  the  source  (located  at  the  origin),  the  pressure  is  an 
initially  spherical  wave  of  the  forn 

p(R,t)  =  (t  >  R/c) 

=  0  (t  <  R/c)  (4) 

ana  R  is  the  range  from  the  source,  while  c  is  the  velocity  of  sound  in  the 
neighborhood  of  the  source. 

For  mathematical  convenience,  we  model  the  source  pulse  [denoted  by 
F(t)j  as  a  linearly  increasing  line  segment  followed  by  a  linearly  decreasing 
line  segment.  Thus,  the  model  source  pulse  is  described  by  the  equations 


F(t)  =  0 

t.  <  0 

(a) 

0  <  t  <  t 
—  —  o 

(b) 

A(t-tx) 

t  <  t  <  t. 
o  —  —  X 

(c) 

=  0 

t  >  t^ 

(d) 

(5) 

and  is  shown  in  Fig.  3* 

Turning  our  attention  bach  to  Eqs.  1  and  2,  we  note  that  the  exponents 
1.13  and  -0.22  appearing  in  Sqs.  1  and  2  would  be  replaced  by  unity  and  zero, 
respectively,  in  an  absorption-free  linear  fluid,  corresponding  to  spherical 
spreading  and  s  time  constant  independent  of  distance."3  The  peculiar  exponents 
appear  due  to  a  combination  of  non-linear  effects  of  the  water  medium  and 
viscous  attenuation.  We  cahe  the  assumption  that  both  these  effects  accumu¬ 


late  with  distance  from  the  source  so  that  tlx*  exponents  1.13  and  -0.22  can 
be  replaced  by  one  and  zero,  respectively,  in  the  neighborhood  of  the  source. 


With  this  assumption,  we  determine  A  by  setting  the  peak  value  of  Eq.  4 


equal  to  the  peak  value  of  (the  modified  form  of)  Eq.  1,  i.e.. 


=  2.25  x  10 


!,»« 


or 


A  =  43  (2.25  X  104)  V1/3 


Setting  the  exponent  -0.22  in  Eq.  2  equal  to  zero,  the  decay  time  in  the 
neighborhood  of  the  source  is  then 

tjj  =  58  (micro-sec) 

or,  for  w  =  0.122  lbs, 

tjj  =  2.87  x  10 secs 

For  convenience,  we  will  henceforth  take  t^  as  10-'*  secs;  it  will  turn  out 
that  the  theoretical  predictions  are  only  weakly  dependent  on  t^. 

Note  from  Eq.  1  that,  for  small  values  of  (t-R/c)/t^,  .•e  have 
(t-R/c) 

T) 

V (t-R/c) 
i7) 


Comparing  the  above  with  Eq,  5c,  we  see  that  if  the  rise  time  of  the 
model  tQ  is  such  that  tQ  «  t^,  as  it  should  be  physically,  Eq.  5c  will 
approximate  the  (early  times)  linear  portion  of  the  decay  of  Eq.  1  if 
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and 


&r5M  « i 


or 

(t-R/c)  «  10~^  secs. 

That  is,  the  model  source  pulse  will  accurately  represent  the  decay  of  (the 
modified  form  of)  Eo.  1  only  for  a  time  interval  up  to  10~^  secs  or  so  after 
occurrence  of  the  peak  value. 

In  the  absence  of  caustics  and  viscosity,  the  pressure  predicted  by 
the  Silbiger  theory  is,  aside  from  a  multiplicative  amplitude  factor  de¬ 
pendent  on  geometry,  a  tine  delayed  replica  of  the  source  pulse. 

After  a  pulse  has  propagated  through  a  caustic,  the  pressure  tine 
2; 

history  is  of  the  fcra 

p(x,y,z,t)  =  X  Pr  f  dT 

where  K  is  a  geometry-dependent  factor,  ?(T)  is  the  time  history  of  the 
source  pulse,  t  is  the  tine  of  observation  and  Pr  denotes  principal  value. 

Proa  the  definition  of  the  principal  value  integral,  it  can  be  shown 
that  the  pressure  field  will  diverge  at  certain  times  if  F(T)  has  any  dis¬ 
continuities.  Therefore,  in  the  absence  of  viscosity,  the  tcodel  source 
pulse  mu5t  have  non-zero  rise  and  decay  times,  lest  the  theory  predict 
infinite  val tes  for  the  pressure  in  the  presence  of  caustics. 
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It  is  shown  in  a  later  section  that  the  effects  of  viscosity  can  be 
incorporated  into  the  theory  by  the  introduction  of  a  complex  propagation 
constant.  Expressed  as  a  Fourier  integral  ever  the  Fourier  spectrum  ?(cu) 


of  the  source  pulse,  the  pressure  appears  as 

p(x,y.z,t)  .  K  f  e-tot-  bKWOW)  ?(o))  ^ 

•OD 

where  cu=2«f  is  the  circular  frequency.  The  presence  of  the  exponential  damp¬ 
ing  factor  can  be  shown  to  insure  a  finite  pressure  at  all  times,  even  in 
the  presence  of  discontinuities  in  F(T),  provided  that  F(T)  is  finite  and  is 
of  finite  duration.  Thus,  if  the  medium  is  assumed  to  be  viscous,  the  model 
source  pulse  may  have  vanishing  rise  and/or  decay  tines. 

As  it  is  well  known  that  sound  waves  propagating  in  water  suffer  viscous 
attenuation,5  two  models  of  the  water  medium  are  considered  for  the  sake  of 
comparison:  one  with  viscosity  included,  the  other  without. 

2.  Modeling  the  Sound  Speed  Profile 

The  sound  speed  profile  is  modeled  by  the  bilinear  profile  shown  in 
Fig.  4a,  along  with  its  associated  ray  diagram  (Fig.  l»b)  corresponding  to  the 
experimental  sound  source  location.  The  model  profile  neglects  the  iso- 
velccity  segment  of  the  experimental  profile  which  originates  at  the  water 
surface,  as  well  as  the  quarry  bottom,  walls  and  water  surface. 

The  caustic  in  the  inhomogeneous  layer  corresponds  to  the  experimental 
caustic.  The  second  caustic  extends  to  infinite  distances  only  because  the 
surface  isovelocity  layer  is  omitted  in  the  model.  If  the  surface  iscvelocity 
layer  were  included,  consideration  of  the  permissible  ray  paths  would  show 
that  the  second  caustic  would  only  extend  a  short  distance  from  the  common 
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point  of  origin  of  the  two  caustics.  Some  ambiguity  is  introduced  by  the 
omission.  On  the  one  hand,  in  the  geometrical  acoustics  limit  the  pressure 
pulse  is  determined  by  sound  which  has  traversed  ray  paths  lying  entirely 
within  the  sloping  portion  of  the  actual  profile,  and  these  critical  ray 
paths  are  reproduced  in  the  bilinear  model.  However,  experimentally,  a  signal 
with  a  distinct  beginning  is  observed  to  precede  the  onset  of  the  caustic- 
focused  pulse.  This  precursor  signal  most  likely  arises  from  sound  which 
has  propagated,  at  least  in  part,  along  the  intersection  between  the  surface 
isovelocity  layer  and  the  sloping  portion  of  the  actual  profile.^  The 
bilinear  model  will  be  unable,  consequently,  to  predict  this  signal  and, 
indeed,  some  uncertainty  in  the  interpretation  of  the  experimental  records 
is  the  result. 

The  omission  of  the  bottom,  walls  and  water  surface  is  justified  since 
a  pulse  reflected  from  any  of  these  surfaces  will  arrive  at  the  observation 
point  much  later  than  the  caustic-focused  pulse,  and  in  practice,  produces 
an  easily  identifiable  trace  in  the  experimental  records. 

The  velocity  profile  of  the  bilinear  model  is  described  mathematically 
by  the  relations 

c(z)  =  cq  (l-az)  z  <  0 

=  c  z  >  0  . 

o 

The  Silbiger  theory  then  contains  three  parameters:  the  gradient  (slope)  of 
the  linearly  increasing  segment  (denoted  by  a),  the  distance  of  the  source 
from  the  origin  (denoted  by  z^)  and  the  sound  speed  at  the  origin  (c  ). 

v*  O 


-8- 


These  parameters  are  fixed  by  requiring  that  the  model  reproduce  (i)  the 
experimental  vertical  distance  between  the  source  and  point  of  observation 
on  the  caustic  {at  the  experimental  horizontal  range),  (ii)  the  experimental 
horizontal  range  (free  the  source)  of  the  origin  of  the  caustics,  and  (iii^ 
the  average  sound  speed  of  the  lower  isovelocity  segment  of  the  experimental 
profile.  A  reasonable  fit  of  the  experimental  situation  was  obtained  with 

z  =  20.8  feet 
o 

a  =  2.707  x  10~3  (ft)'1 
cq  -  U715  ft/sec 

Smaller  values  of  the  parameter  a  yielded  significant  disagreement  with  the 
experimental  range  to  the  caustic  origin;  larger  values  produce  e  large 
discrepancy  with  the  experimental  vertical  distance  between  source  and  caustic. 
A  comparison  between  the  mod^l  and  the  experiment  is  shown  in  Table  1,  along 
with  some  other  choices  of  model  parameters  to  provide  a  comparison.  The 
quantity  (zq-zc)  is  the  vertical  distance  between  the  source  and  the  caustic. 

R  is  the  horizontal  range. 

Note  that  the  horizontal  range  to  the  origin  of  the  caustic  is  not  veil 

defined  experimentally  (see  Table  1  and  Fig.  2).  According  to  the  bilinear 

model,  ZQ~ZC  should  decrease  nonotonically  as  a  function  of  horizontal  range, 

whereas  in  Table  1  the  experimental  value  cf  z  -z  is  non-monotonic  for  the 

o  c 

closest  ranges.  Our  choice  of  model  parameters  yields  a  range  of  227.7  feet. 
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•  ir.-l 


2  =  20.8  ft  2  =  20.0  ft  zn  «  18.0  ft  j  z_  =  17.0  ft 


8=2.707*10“3ft“1  a=2 .  TOTrlO’^ft"!1  a=2c02lacl0“3ft‘1  a=1.3i»lxlO'3ft_1 


z  -z„ 
o  c 


z -z„ 
0  c 


2  “Z„ 

o  c 


Vzc 


193.1!  28.6 
209.8  27.8 
210.1  28.0 

215.3  27.1 


I 

227.7 

23.3 

232.2  25.7 

231.7 

23.2 

21*3.5  2^.9 

21*3-7 

l 

22.8 

257.6  21*. 1 

j  257.7 

22.5 

I 

275.3  23.1 

1 

!  275.7 

22.2 

293.0  22.1 

1 

297.7 

21.9 

331.9  20.0 

;  332.0 

21.7 

378.1  17.6 

1  ■ 

378.0 

21.1* 

389.9  16.7 

3«3.0  j 

21.1* 

l*l6.p  15.0 

!  1*00.0 

21.1* 

2V».3  |  21.8  2lt5.0  20.2 

259.3  '  21.5  258.0 


258.0  j  19.8 
19.5 


299.3  :  21.0 


3«  Pressure  Tine  History  on  Caustic:  Inviscid  Medium  Model 

For  points  on  the  caustic  (in  a  non-vjscous  medium),  the  time  history 

7 

of  the  pressure  is  given  by 


p(vvt) 


8k 


c(o) 


■1/3 


rAj-^F  tn^O-S?]2 


J  $(<»)  exp(lo)[Wc/c(o)-t-  |/|cuf3)  &i> 


(6) 


where  r  and  z are  the  cylindrical  coordinates  of  the  point  of  observation 
on  the  caustic,  c(o)  is  the  sound  speed  at  the  origin,  §  is  the  sine  of  the 
angle  with  respect  to  the  vertical  at  which  the  ray  touching  the  caustic  et 
the  observation  point  leaves  the  source,  n(zc) [n(zQ)] is  the  index  of  refrac¬ 
tion  at  the  caustic  (source),  w"  is  the  third  derivative  with  respect  to  % 
of  the  phase  function  V(§,z,r)  evaluated  at  'c>zc>rc>  vhere 


W(S,z,r)  =  §r  +  j*  c  ;'n2(z)-?2  dz  +  J  /n2( z)~^  dz 
z_  '  z„  * 


and  z_  is  the  turning  point  for  the  ray  characterized  by  §,  Ai(o)  is  the 
Airy  function  evaluated  at  zero  argument,  r*(co)  is  the  Fourier  transform  of 
F(t),  defined  by  the  relation 


+50 


¥(<d)  -  J  F(t)  e^  dt 
and  for  the  model  source  pulse  (Fig.  3)  is  given  by 


~  4  'e^o  -1)  inti  iutQ 

VU)  =  JL  rve  .  9.  rJJ.  +  « ±  .  « — — 1 

{  )  2  1  fc-  Vh.  VV 


CD 
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+ 


(7) 


to  is  the  circular  frequency,  is  W(§c,z^,rc)  and  t  is  the  tins.  Using 
Eq.  7,  the  pressure  can  be  written  in  closed  form: 


=  c  <?  Klrl/3 


2jrS„  c(o) 


-1/3 


rctn2(zc)-rcBl-§2cP 


Ai(o) 


*1 


W 

t - 7— r  -  t 

c(o)  O 


5/6 


(t‘  ‘  V<  ^+1)}  +  t" 


t- 


{( /T-l  +  SON 

5/6 


{f/^-l)+SGN(t-  ( V3  4-1)} 


-  XrM;-  £7  -tJtfjH)} 


(8) 


where  SGN(x)  =  1  x  >  0 

=  -1  x  <  0 

The  experimental  data  of  Barash  is  presented  not  in  terms  of  the 
pressure  directly,  but  rather  in  terms  of  the  ratio  of  the  observed  pressure 
pulse  to  the  peak,  value  of  the  pressure  pulse  which  would  be  observed  if  the 
water  had  a  uniform  velocity.  This  ratio  is  called  the  amplification 
factor.  For  isovelocity  water,  the  peak  value  of  the  pressure  pulse  would 


where  w  is  the  weight  (in  pounds)  of  the  explosive  and  r  ,  zc  and  z^  are 
measured  in  feet. 
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Denoting  the  last  bracketed  expression  on  the  right  hand  side  of  2q.  6 


by  f(t,t^,tQ),  the  amplification  factor  can  be  written  as 


P(r  .z  t) 


■f-2 —  s  *  f(t,tx,t  ) 

iso 


(9) 


The  factor  FPP  is  determined  by  the  location  of  the  point  of  observation 
and  the  parameters  of  the  bilinear  model.  FPP,  as  a  function  of  the  model 
parameters,  is  given  in  Table  2.  The  horizontal  range  (Rc)  to  the  caustic 
origin  is  also  shown.  FPP  is  strongly  dependent  on  the  slope  (a),  r 

c 

vari  '  with  both  z  and  a.  As  mentioned  before,  the  first  set  of  parameters 
o 

offers  the  most  reasonable  fit  to  the  experimental  caustic  curve. 

The  maximum  value  of  the  function  f(t,t^,t  )  is  shown  in  Table  3  along 
with  the  corresponding  maximum  value  of  the  amplification  factor,  for  a 
series  of  values  of  t^  and  tQ.  The  observed  maximum  amplification  factor 
is  b.95.  All  the  predicted  values  are  too  high. 

If  the  precursor  pressure  and  the  caustic-focused  pressure  combine 
additively  to  yield  the  experimental  record,  the  experimental  amplification 
factor  would  be  reduced  to  about  3*3-  It  is  not  clear,  however,  if  this 
interpretation  is  correct  but  a  theoretical  study  cf  precursors  is  in 
progress.  In  any  case,  the  predicted  values  will  still  be  too  high. 

The  function  f(t,t^,tQ)  is  shown  in  Fig.  5.  Values  of  f(t,t^,t  )  for 
ties  intervals  larger  than  about  10  ^  seconds  on  either  side  of  the  origin 
are  not  expected  to  be  accurate  because:  (l)  the  model  source  pulse  is  valid 
only  up  to  times  of  10-^  seconds  or  so,  and  (2)  the  high-frequency  approxi¬ 
mations  made  in  the  theory  itself.  These  latter  limitations  have  not  been 


Table  3 

MAXIMUM  VALUES  OF  AND  THE  AMPLITUDE  FACTOR 

AS  HJNCTIOIB  OF  t,  AND  t 
1  o 

t.(sec)  tQ(sec)  KAXXMUM  AMP  FACTOR 


io'5 

io'7 

37*37 

9-6 

IO-5 

10-8 

61.22 

15.8 

3.0“5 

10"S 

95*94 

24.7 

10-“ 

IO'6 

25-4 

6.5 

IO-4 

10-8 

69.2 

17.? 

10'4 

2.1  X  10~1T 

2091 

539 

10-“ 

3  X  IO-18 

3341 

660 
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investigated  thoroughly  enough  to  justify  a  discussion  of  theca  here.  Hovever, 
a  study  of  the  limitations  of  the  caustic  boundary  layer  formulae  is  in 
progress.  The  slowly  rising  signal  occurring  prior  to  t=0  has  no  definite 
time  of  onset  and  may  be  associated  with  the  high  frequency  nature  of  some 
of  the  approximations  employed  in  the  theory.  Its  significance  has  yet  to 
be  studied  in  detail. 

The  experimentally  observed  time  history  (of  the  amplification  factor) 

is  shown  in  Fig.  6.  The  precursor  signal  precedes  the  actual  focused  pulse. 

The  theory,  as  presently  formulated,  cannot  predict  this  signal. 

The  details  of  the  experimental  record  Eire,  to  some  extent,  unreliable 

in  the  s^nse  that  a  sudden  rise  in  pressure  was  recorded  with  a  10  micro- 
2 

second  rise  time.  Hence,  although  it  appears  in  Fig.  6  as  though  the  ex¬ 
perimentally  observed  rise  and  decay  times  of  the  focused  pulse  are  both  on 
the  order  of  10  ^  secs,  they  cay  in  fact  be  shorter  than  that. 

At  any  rate,  we  can  conclude  from  Table  3  that  the  rise  tine  required 
by  the  theory  to  reproduce  the  experimental  peak  amplification  factor  is 
unrealistically  large,  since  it  must  be  longer  than  10  secs,  for 
t,  =  10  5  secs. 

n.  Pressure  Tims  History  at  Caustic:  Viscous  Medium  Model 

The  tin?  dependent  wave  equation,  including  Viscosity,  is  given  by^ 

•  A  >2 

v  p(x,y,z,t)  +  ^2p(x,y,z  t)  — - x  p(x,y,2.t) 

c2(z)  c(z)  dt2 

=  -  6(x)5(y)S(z)  [F(t)  +  §-  F(t)3 

c  (o) 


where  v  is  the  kineoatic  viscosity  coefficient  defined  as 

k 

n  *  ?  U 


tj  =  coefficient  of  bulk  viscosity 
p  =  coefficient  of  shear  viscosity 
p  =  density 

For  a  sinusoidal  source  (e  time  dependence)  of  strength  ?{co),  the 
viscous  wave  equation  reduces  to 


2  2 

+  "'"'dffi  '  P  *  *  ^(®)  6(x)6(y)f(z) 

C2(z) 


and  the  variation  of  c(z)  in  the  spatial  dona  in  of  interest  is  snarl  enough, 

then  in  the  factor  [1  -  — x ,  c(z)  cr>y  be  replaced  by  e(o)  and  an  effec- 

c  (-.) 

tive  propagation  constant  defined  as  follows 


c2(o) 


Introducing  the  Hankel  transfora  pair 


f(§,z)  =  5  P(r,z)  J (k §r)  r  dr 
0 
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(13) 


p(r,z)  =  J  f(§,z)  JQ(keff,  5r)  k2ff  §  d§ 


the  transform  of  Eq.  10  becomes 


dz 


f(§,z)  +  k “ff  [n“(z)  -  l2)  f(S,z)  =  - 


6(z) 


(1*) 


Proceeding  free  Sas.  12  and  13,  the  development  of  the  caustic  boundary 
layer  formalism  proceeds  as  in  the  non-visccus  case,  the  only  difference 
being  the  appearance  of  koff>  instead  of  k. 

The  pressure  time  his^vry  on  the  caustic  is  then  given  by 

p(re,zc,t)  =  f  c/ c{o)-t}«*  l/G 


where 


w  =h 


2a  §c  c(c) 


+  complex  conjugate] 

-1/3 


rc(n2(Ic)-^)^(»2Cz0)-^)^ 


(IKI)'1'3 


and 


a  = 


2cJ(o) 


For  the  model  source  pulse,  the  pressure  time  history  is  given  by 

_/_  _  ^  (AA)-Ai(o)*A.V2r  *1  *c  _  ,  1  *c 

p(rc- — {ftvTt  H[^y  ‘  mo3-  r  ' t3 

0X0  o 


-  TJpry  Hi^  ■  wt!)J 


(15) 
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where 


H(t)  =  -  |  S[3'7/l2r(7/12)  ^ 

♦  «-«/“  A  # 

+  |  »C3-1/12  r(l/12)  A  4;  §,•  -  J|£> 

-t  ^r(,/12)A(y;-|)! 


and 

t?  W 

c  =  c 
2c3(o) 

Using  the  asymptotic  form  of  -F^  valid  for  t  lc/lj3  »  1,  it  can  be 
verified  that:  (a)  Eq.  15  reduces  tc  the  non-viscous  pressure  in  the  limit 
■„*-*  o  (no  viscosity)  and  (b)  Eq.  15  reduces  to  the  non-viscous  pressure  in 
the  limit  to  |  tQ j2/i*8  »  1,  .tjJ^/4c  >:>  Toe  first  limit  provides  a  check 
on  Eq.  15.  The  second  limit  implies  that  for  sufficiently  long  rise  and 
decay  times,  viscous  attenuation  has  no  effect  on  tie  pressure  time  history. 
This  is  reasonable  since  long  rise  and  decay  times  means  there  is  very  little 
high  frequency  content  in  the  Fourier  spectrum  of  F(t).  This  being  the  case, 
it  is  Immaterial  whether  the  high  frequency  components  get  damped  out  or  not. 

*5 

(Recall  that  viscous  attenuation  goes  as  e  con3tant  “  ^ 

A  function  analogous  to  f(t,t^,tQ)  of  Eq.  9  can  be  defined  for  the 

viscous  model.  For  a  model  source  pulse  with  zero  rise  time  and  tj=10  ^ 

seconds .  the  function  [denoted  by  fv(t)]  is  shown  in  Fig,  7.  The  viscosity 
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related  parameters  are  for  fresh  water v  u  =  0.01  g/cm-sec  and  tj=2.81  m* 
c(o)  was  tahen  as  1.5  x  10^  cm/sec  ana  as  300  ft,  the  distance  along  the 
ray  path  to  the  point  of  observation.  The  amplification  is  then  taxen  as 
f^(t)  tires  FFP  as  defined  in  Eq.  9* 

As  in  the  non-viscous  case,  the  computed  tire  history  has  a  limited  time 

domain  of  validity.  A  numerical  evaluation  of  the  inequality  of  Eq.  11  shews 

2  10 
that  cou/c  (z)  is  on  the  order  of  0.1  for  a  frequency  of  10  Hz.  Therefore 

the  computed  time  history  should  be  unreliable  for  times  shorter  than  about 
10  seconds.  As  mentioned  before  in  connection  with  the  inviscid  medium, 
the  model  source  pulse  becomes  unreliable  for  times  longer  than  10-^  seconds 
or  so.  Therefore,  the  computed  time  history  should  become  erroneous  for  tines 
(beyond  the  origin)  of  the  same  magnitude.  In  addition,  the  theory  is  based 
on  a  high  frequency  approximation  which  will  render  the  time  history  reliable 
for  only  a  short  time  interval  after  onset.  However,  a  quantitative  descrip¬ 
tion  of  this  source  of  error  is  not  yet  available. 

For  the  model  source  pulse  (Fig.  3)  with  tQ=0,t^=10  5  secs,  the  theoreti¬ 
cal  peafc.  amplification  factor  is  6.1,  as  compared  to  the  experimental  value 

of  4.95* 

If  the  precursor  pressure  and  the  caustic-focused  pressure  oabined 
additively  to  yield  the  experimental  record,  then  the  experimental  amplifi¬ 
cation  factor  would  be  about  3.3. 

In  either  case,  the  predicted  value  is  too  high  by  20  to  100  percent.  If 
a  non- zero  rise  tire  in  the  model  source  pulse  were  assumed,  the  theoretical 
prediction  should  decrease.  It  is  possible  that  a  more  reasonable  prediction 
could  be  made  with  a  non-zero  yet  physically  reasonable  rise  time,  but  this 
has  non  oeen  attempted. 
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HI.  CCKCUJSI0S3 

Assailing  a  source  with  finite  rise  time  and  an  inviscid  medium,  the 
computed  e nullification  factors  are  much  too  high  if  a  physically  acceptable 
rise  time  is  assumed. 

For  a  source  with  zero  rise  time  radiating  into  a  viscous  medium,  the 
theoretically  computed  nullification  factor  is  20  percent  to  100  percent  too 
high.  A  non-zero  rise  time  might  reduce  the  discrepancy. 
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Fig.  1.  Experimental  sound  velocity  profile.  (Courtesy  of  R.  H.  Barash  and  R,  Tnran, 
U.S.  Naval  Ordnance  Laboratory.) 


INHOMOGENEOUS  LAYER 


i’ J-ji? •  *' •  Ma  bhoinatical  prediction  of  pressure  pulse  shape 
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